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£SJ , Abstract 

O j We consider two limiting regimes, the large-spin and the mean-field limit, for the dynamical 

evolution of quantum spin systems. We prove that, in these limits, the time evolution of a class of 
quantum spin systems is determined by a corresponding Hamiltonian dynamics of classical spins. 
qq , This result can be viewed as a Egorov-type theorem. We extend our results to the thermodynamic 

■ limit of lattice spin systems and continuum domains of infinite size, and we study the time 
evolution of coherent spin states in these limiting regimes. 

1 Introduction 

The purpose of this letter is to study classical limits of quantum spin systems. Work in this direction 
was undertaken already at the beginning of the seventies. Yet, most of the mathematical results 
established so far only concern time-independent aspects, such as the classical limit of quantum 
partition functions for spin systems [H [5] . Here we consider the dynamical evolution of quantum spin 
systems in limiting regimes; see also [7]. In particular, we discuss i) the large-spin limit and ii) the 
mean- field /continuum limit. As our main results, we prove that the time evolution of a large class of 
quantum spin systems approaches the time evolution of classical spins. Our results can be regarded 
as Egorov-type theorems, asserting that quantization commutes with time evolution in the classical 
■^J" ' limit; see [T] for a similar result on classical and quantum Bose gases. Along the way, we discuss 

, thermodynamic limits and the time evolution of coherent spin states, in the two limits mentioned 

■ above. 

£"■■» , An example of an evolution equation for classical spins is the Landau- Lifshitz equation 



O 



dtM = M A ffcx(M) , (1) 



which is widely used in the study of ferromagnetism. Here M = M(t, x) G S 2 denotes a classical spin 
field with values on the unit sphere, a 
the exchange field is H CX (M) = J AM 
constant. Equation ([I]) then becomes 



field with values on the unit sphere, and A stands for the vector product in K 3 . A standard choice for 
the exchange field is H CX (M) = J AM, where A denotes the Laplacian and J is the exchange coupling 



d t M = JM A AM . (2) 

This form of the Landau-Lifshitz equation has been studied in the mathematical literature; see for 
instance [5J [3] and references given there. In physical terms, @ describes the dynamics of spin waves 
in a ferromagnet with nearest neighbor exchange interactions in a classical regime; see [B]. 
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In this paper we consider the Landau-Lifshitz equation with an exchange field H ex (M) given by 
an integral operator applied to M, and generalizations thereof. Equation (|TJ) then takes the form 

d t M(t, x) = M (t, x) A J J{x, y) M (t, y) dy . (3) 

The integral kernel J(x, y) describes the exchange interactions between classical spins beyond the 
nearest-neighbor approximation in the continuum limit. A formal argument on how to derive ^ from 
is given in a remark in Sect. 13.41 

Our paper is organized as follows. In Sect. O we study the dynamics of finite lattice systems of 
quantum spins in the limit where their spin s approaches oo. The main result of Sect. [2] is formulated 
in Theorem [T] below. In order to prepare the ground for this theorem and its proof, we first introduce 
a class of Hamilton functions for classical spins and define their quantization by means of a normal- 
ordering prescription. At the end of Sect. [2 we pass to the thermodynamic limit, and we discuss the 
time evolution of coherent spin states. 

In Sect. [3l we present a similar analysis for the mean-field limit of quantum spin systems defined 
on a lattice with spacing ft, > in the continuum limit, h — > 0. The main result of Sect. [3] is stated in 
Theorem H] below. 



2 Large- Spin Limit 

2.1 A system of classical spins 

Let A be a finite subset of the lattice Z d (or any other lattice). A classical spin system on A is 
described in terms of the finite-dimensional phase space 

r A := n §2 > 

i. e. we associate an element M(x) of the unit two-sphere § 2 c M 3 with each site x S A. The phase 
space Ta is conveniently coordinatized as follows. For each site x € A, let (Mi(x), M^x), M^{x)) 
denote the three cartesian components of a unit vector M{x) <E S 2 , and define the complex coordinate 
functions (M+(x), M z (x), M_(x)) on S 2 by 

. . M 1 (x)±iM 2 (x) „, , , „, , . 

We define a Poisson structur^ on Ta by setting 

{Mi(x),Mj(y)} = is ijk S(x,y)M k (x), (4) 

Here 5(x,y) stands for the Kronecker delta, and the indices i,j, k run through the index set / := 
}, where the symbol e^-fc is defined as e± TZ — ±1, e± z ± — e z ±± = ±1, and iijk = 
otherwise. 

For our purposes it is convenient (but not necessary) to replace § 2 with the closed unit ball 
Bi(0) C M 3 . To this end, we introduce a larger "phase space" (a Poisson manifold) 

xeA 

equipped with the Z°°-norm. The algebra 

yS A :=C[{Mi(x) : i e I, x e A}] 

of complex polynomials is a Poisson algebra with Poisson bracket determined by We equip *Pa 
with the norm ||-A||oo := sup Mg5 ; A |A(M)|, and we denote its norm closure by 21a- Note that, by the 
Stone- Weierstrass theorem, 21a is the algebra of continuous complex- valued functions on Sa ■ 

1 Actually Ta is symplectic, with symplectic structure determined by the usual one on S 2 . 
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A fairly general class of Hamilton functions H\ on Ha may be described as follows: We asso- 
ciate with each multi-index a £ N /xZ satisfying |a| := J2 x ez d YliPi a ii x ) < °o a complex number 
V(a). Using the trivial embedding N IxA C N /xZd defined by adjoining zeroes, we consider Hamilton 
functions of the form 

H A ■= v ( a ) Ma - (5) 

In order to obtain a real-valued Ha, we require that V(a) — V(a), where the "conjugation" a of a 
multi-index a is defined as cii(x) :— aj(x), with ~ : (+,z,— ) i— > (— ,z,+). Furthermore, we impose 
the following bound on the interaction potential^ 

|| V|| := ]T SU P E Hx)\\V(a)\e n < oo . (6) 

\a\—n 

It is then easy to see that the series converges in norm and that the set of allowed interaction 
potentials V is a Banach space. The Hamiltonian equation of motion reads A — {Ha, A}, for any 
observable A £ 21a- In particular, a straightforward calculation yields 

j t M t (t, x ) = J2 V{a)^ii Jikaj {x)M°>- s 'V+ s >>M{t), (7) 

a£S ,xJ j,k 

where the multi-index 5i(x) is defined by [5i(x)]j(y) '■= 5ij5(x,y). 

We record the following well-posedness result for the dynamics generated by the class of Hamilto- 
nians introduced above. 

Lemma 1. Let A be a (possibly infinite) subset o{1 d . Let Mq £ Ha- Then the Hamiltonian equation 
has a unique global-in-time solution M £ C 1 (R, Sa) that satisfies M(0) = Mq. Moreover, the 
solution M depends continuously on the initial condition Mq, and we have the pointwise conservation 
law \M(t,x)\ = |M(0,a:)| for allt. 

Proof. Local-in-time existence and uniqueness follows from a simple contraction mapping argument 
for the integral equation associated with (|7|). We omit the details. Also, continuous dependence on 
Mo follows from standard arguments. Finally, the claim that |M(0,x)| = \M(t, x)\ for all t can be 
easily verified by using ([7]), which implies that 4fM(t,x) is perpendicular to M(t,x). □ 

Remarks. 1. In what follows, we denote the flow map Mo i— > M[t) by 4> A . Note that, under our 
assumptions, ([7]) also makes sense for infinite A C Z d , whereas the Hamiltonian Ha does not have a 
limit when |A| — ► oo. 

2. The last statement implies that the magnitude of each spin remains constant in time, i.e. the 
spins precess. In particular, if Mo £ Ta, it follows that M(t) £ Ta for all t. Mathematically, this is 
simply the statement that the symplectic leaves of the Poisson manifold Ea remain invariant under 
the Hamiltonian flow. 

3. Time-dependent potentials V(t, a) may be treated without additional complications, provided that 
the map t i— » V(t) is continuous (in the above norm) and sup^. in © is replaced by sup^j. The 
weaker assumption that 1 1— > V(t, a) is continuous for all a implies Lemma [T] with the slightly weaker 
statement that M £ C(K, Ha) is a classical solution of (J7J). 

Example. Consider the Hamiltonian 

where M(x) = (Mi(x), M2(x), M^{x)). Here h(t,x) £ M 3 is an "external magnetic field" satisfying 
sup tgR X £.z,d\h(t, x)\ < oo. We also require the map t t— * h(t,x) to be continuous for all x £ 7L d . 

2 Note that this condition may be weakened by replacing e™ with e rn , for any r > 0. It may be checked that this 
does not affect the following results. 
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The exchange coupling J : Z d x Z d — > R is assumed to be symmetric and to satisfy J(x, x) = for 
all x. Finally we assume, in accordance with condition ([H]), that sup^g^ *^2 ye id\J(x, y)\ < oo. The 
corresponding equation of motion for M(t, x) is given by 



dt 



M(t,x) = M(t,x) A 



h(t, x ) + J2 J ( x >y) M ( t >y) 



(9) 



i. e. the Landau- Lifschitz equation for a classical lattice spin system. 



2.2 A system of quantum spins 

In this section we formulate the quantum analogue of the system of classical spins from the previous 
section. We associate with each point x £ 1 d a finite-dimensional Hilbert space J^f = := C 2s+1 
describing a quantum-mechanical spin of magnitude s. (Here, and in the following, we refrain from 
displaying the explicit s-dependence whenever it is not needed.) Furthermore, we associate with each 
finite set Actf the product space J?a &> x eA<^x> an< ^ we define the algebra 21a as the algebra of 
(bounded) operators on J4?\, equipped with the operator norm || • ||. 

The spins are represented on Jff\ by a family {Si(x) : i = 1,2, 3, x G A} of operators, where 
Si(x) is the i'th generator of the spin-s-representation of su(2) on Jf? x , rescaled by s _1 . In analogy to 
the complex coordinatization of the classical phase space Ta in the previous section, we replace the 
operators (Si(x), S2(x), S3 (a:)) with (S+(x), S z (x), S-(x)) as follows: 

S ± ( x ) := gjM±iW t g z (x) := S 3 (x) , for all x e Z d . 
v2 

An easy calculation yields |p±(a;)| < 1 and ([^(a;)! = 1 if s > 1, and ||5±(x)|| — \/2 and ||5 z (x)|| = 1 
if s = 1/2. Furthermore one finds the fundamental commutation relations 

[Si(x),Sj(y)] = -iij k 5{x,y)Sk{x) (10) 

with i, j, k S /. 



2.3 Quantization 

In order to quantize polynomials in ^Pa we need a concept of normal ordering. We say that a monomial 
5*11 {x\) ■ ■ ■ Si p (x p ) is normal- ordered if ik < h => k < I, where < is defined on / through + < z < — . 
We then define normal-ordering by 

: (xi) ■ ■ ■ S ip (x p ) : = S ia{1) (x (T (i)) ■ • • S itr(p) (x CT ( p )) , 

where a € S p is a permutation such that the monomial on the right side is normal-ordered. Next, we 
define quantization ~ : ^\ 21a by setting 

(M^xi) ■ ■ ■ M ip (x p )) = : S il (xi) ■ ■ ■ S lp (x p ) : 

and by linearity of ^ . We set 1 = 1. Note that, by definition, ^ is a linear map (but, of course, not 
an algebra homomorphism) and satisfies {A)* = A. 

2.4 Dynamics in the large-spin limit 

For each finite A C Z d we define the Hamiltonian H\ as the quantization of Ha- More precisely, we 
quantize ([5|) term by term and note that the resulting series converges in operator norm. Because 
is real, the operator H\ is self-adjoint on the finite-dimensional Hilbert space Jff\ and generates a 
one-parameter group of unitary propagators U s (t;H\) (equal to e ~ lsHAt if H\ is time-independent). 
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We introduce the short-hand notation 

a A A := A o <p A , A G 21a > 

aiA := U s (t;H A )*AU s (t;H A ), Ae&a, 

where <j> K is the Hamiltonian flow on 5a. Note that both a A and a\ are norm-preserving. 

We are now able to state and prove our main result for the case of a finite lattice A. Roughly it 
states that time evolution and quantization commute in the s — > oo limit. This is a Egorov-type result. 



Theorem 1. Let A e <}3 A and e > 0. Then there exists a function A(t) £ ^a such that 



sup||a A A- ^(iJlloo < e, 
teR 



and, for any t€l, 



|a A A-A(f)|| < e 



C(e,t,A) 



(11) 



(12) 



where C(e,t,A) is independent of A. 



Proof. Without loss of generality we assume that A = M@ for some e N /xA . For simplicity of 
notation we also assume, here and in the following proofs, that Ha is time-independent. Consider the 
Lie-Schwinger series for the time evolution of the classical spin system, 



00 A 



(13) 



where {H A ,Af = {h a ,{Ha,A} (1 X) } and {Ha,A} (0) = A. In order to compute the nested 
Poisson brackets we observe that 



xeA i.j.kei 

as can be seen after a short calculation. Iterating this identity yields 

{ha,a } ^ = i« £ E E E E 

a 1 , .••>«' Xi,...,Xi fei,...,fei 
' I r 9-1 

11^^ V ( a9 ) a l( X l) /J + EK" ~ *Jr(*r) +*fcr(a:r)) 

,g— 1 L r— 1 



(14) 



M /3+T, l r=1 (a r -S ir (x r )S jr (x r )+S kr (x r )) _ ^ 

In order to estimate this series, we recall that HM 7 ^ < 1 and rewrite it by using that 

oo 

E - E E ••• E • 

a 1 , ...,a l ni,...ni=l \ct 1 \=ni |a ! |=n; 

We then proceed recursively, starting with the sum over a l ,xi,ii,ji,ki and, at each step, using that 



E EEi^i^^)^^)!^^)! ^ 

|a|=n x i,j,k 



(n) 



where 



(") := sup E \ V ( a )\ K x ) 
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In this manner we find that 

(Oil 



{Ha, A}' 



* E 

oo 

< E (1/ 

oo 

E 

ni,...,m=i 
= l!eW||V||'. 



+ n 1 )---(\f3\+n 1 + ---+n l - 1 )\\V\\ ( - n ^ 

+ ni + --- + ni)'ll^ll (ni) ---H^II (ni) 
|/3|+niH h"i ||\Z||( ni ) . . . ||V||(" ! ) 



Ml 



(ni) 



Thus, for |i| < H^ll- 1 , the series (j!3|) converges in norm, and an analogous estimate of the remainder 
of the Lie-Schwinger expansion of a A A shows that l|13p equals a A A. As all estimates are independent 
of A, the convergence is uniform in A. 

The quantum-mechanical case is similar. Consider the Lie-Schwinger series for the time evolution 
of the quantum spin system: 

oo 



(16) 



where [.Hajj?]^ = H\, [H A ,A\^ 1 ^1 and [iZA,^4]^ = A. In order to estimate the multiple com- 
mutators, we remark that, from Q and (|10[) and since both { • , • } and is[- , • ] are derivations in 
both arguments, we see that (is) l \H\,A\ is equal to the expression obtained from \H\,A\ by 
reordering the terms appropriately and by replacing Mi(x) with Si(x). In particular (assuming s > 1) 

(is) l [H A ,A\ il) \\ < lle?\\V\\ l , 

and we deduce exactly as above that (fT6|) equals a A A for t < ||V|| _1 . 

To show the claim of the theorem for |i| < ||V|| _1 we first remark that a A (A) is well-defined 
through its convergent power series expansion. Now as shown above, each term of (is) 1 \H A , A\ , as a 
polynomial in 21a, is equal to a reordering of the corresponding term of {H A , A\ ' . If P is a monomial 
(with coefficient 1) of degree p in the generating variables {Si(x)} and P a monomial obtained from 
P by any reordering of terms, the commutation relations (|10[) imply that 



IP-PII < 



Thus 



({H A ,A} {1) ) = (isYlH^Al^+Ru 



where, recalling the expression (fT5| and the estimates following it, we see that the "loop terms" Ri 
are bounded by 



\Ri\\ < 



1 



E 



+ n 1 + --- + n ; )' +2 ||^|| ( " l) ---||T/|| (,li) 



Tlx , ...,ni = l 



< 



(l + 2)\ 



V\\ l 



Therefore, if \t\ < \\V\\ 



\a A A-^A\\ < e ^ljr { i + 2)(i + m\\v\\y < ^1 

1=0 



G 



where C(t, A) is independent of A. 

In order to extend the result to arbitrary times we proceed by iteration. The crucial observations 
that enable this process are that the convergence radius is independent of and a A ,a A are 

norm-preserving. Let t G K and choose v G N such that r := tjv satisfies |r| < In order to 

iterate we need to introduce a cutoff in the series (Q~3|) and (fT5|) . The series (fT3|) consists of an infinite 
sum of terms in <Pa which are be indexed by (7, a 1 , . . . , a 1 , X\, . . . , xi). Now let e > be given. Since 
the series converges in norm there is a finite subset 

oo 

B 1 = B x {e) C {(l,a 1 ,...,a l ,x u ...,x l )} = \J(N IxA ) 1 x A 1 

such that the norm of the series restricted to the complement of B\ is smaller than e/v. This 
induces a splitting a A A = a Bl A + a B ^A (in self-explanatory notation), such that a Bl A G *}3a and 

lla^Aljoo < e/v. Similarly, one splits a A A = a Bl A + a B ^A where, after an eventual increase of B\, 

\\a Bi A\\ < e/v. 

Now we use the above result for |t| < \\V\\ 1 : 

alA = a Bl A + — + a B c A 

s 

where R\ is some bounded operator. Since a Bl A G *Pa we may repeat the process on the time interval 
[r,2r]: 



^ g;T ^ 

a A a A A = a A a B A-\ \-a A a B <=A 

= (as 2 a Bl A) + a Bl A + a T A a B ^A 



i? 2 + «a Ri 



Continuing in this manner one sees that, since a A and a A are norm-preserving, A(t) := a Bu ■ ■ ■ a Bl A G 
<Pa satisfies 

WaiA-Ait)^ < e, 



as well as 

112^-^)11 < £ + 

□ 



2.5 The thermodynamic limit 

The above analysis was done for a finite subset A, but the observed uniformity in A allows for a 
statement of the result directly in limit A = Z d . We pause to describe how this works. 

Concentrate first on the quantum case. If Ai C A2, an operator Ai G 21a 1 may be identified in 
the usual fashion with an operator A2 G 21a 2 by setting A2 = A\ ® 1a 2 \Ai • We shall tacitly make 
use of this identification in the following. It induces the norm-preserving mapping 21ai - > 21a 2 °f the 
abstract C*-algebras and the isotony relation 21aj C 21a 2 ■ Observables of the quantum spin system 
in the thermodynamic limit are then elements of the quasi-local algebra 

a := \/ 

AC1 d finite 

which is the C*-algebra defined as the closure of the normed algebra generated by the union of all 
SIa's, where A is finite. The spins are represented on 21 by a family {Si(x) : i G /, x G Z y } of 
operators. 

The dynamics of the system is determined by a one-parameter group (S')t 6 R of automorphisms of 
21. Its existence is a corollary of the proof of Theorem [TJ 
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Lemma 2. Let A £ 21a for some finite Aq C Z d and t £ 1. Then the following limit exists in the 
norm sense: 

lim a l K A =: a* A, 

A^oo 

where A — > oo means that A eventually contains every finite subset. By continuity this extends to a 
strongly continuous one-parameter group (a )tgg of automorphisms o/2l. 

Proof. For |t| < the series (fTBj) is bounded in norm, uniformly in A, so to show convergence of 

the series it suffices to show the convergence of \H\, A\ ^ for each I £ N, which is an easy exercise. 

Thus a* A is well-defined for any polynomial A. By continuity, a* extends to an automorphism 
of 21. Since a 1 A £ 21 and a* is a one-parameter group, we may extend it to all times by iteration. 
Strong continuity follows since S*^4, for small t and polynomial A, is defined through a convergent 
power series: 

lim HSU- -411 = 0. 

By continuity, this remains true for all A £ 21. □ 

For classical spin systems we recall that, for finite A, we have 21a = CfQzeA -Bi(O); C), a C*- 
algebra under || • H^. As above, for Ai C A2, we identify A\ £ 21ai with a function A2 £ 21a 2 by 
setting A 2 — A\ ® 1a 2 \Ai- We thus get a norm-preserving mapping 21a! — > 21a 2 01 the abstract 
C*-algebras and the relation 21a 1 C 21a 2 ■ Define the classical quasi-local algebra as 

21 := \/ 2l A • 

AcZ d finite 

Note that 21 is equal to the space of continuous complex functions on Y\ xeZ d -Bi(O), equipped with the 
product topology (this is an immediate consequence of the Tychonoff and Stone- Weierstrass theorems). 

The spins are represented on 21 by a family {Mi(x) : i € /, x E W} of functions. Existence of 
the dynamics follows exactly as above. 

Lemma 3. Let A € 21a for some finite Ao C 1 d and t € R. Then the following limit exists in || • W^: 

lim a' A A =: a* A. 

A^oo 

By continuity this extends to a strongly continuous one-parameter group (a')tgM of automorphisms of 
21. Furthermore, a 1 A — A o <f> 1 } where <f> 1 — (/>| d is the Landau- Lifschitz flow defined in the remark 
after Lemma\]\ 

Now set <p := C[{Mi(x) : i€l,x€ Z d }}. Then the proof of Theorem □ yields the following 
Theorem 2. Let A £ <p and e > 0. Then there exists a function A(t) £ *p such that 

supHaM- A(t)||oo < £, (17) 
teR 

and, for any i£E, 

WafA-A&W < £+ C M' A ) . (1 8) 
n n s 

Remark. In particular, the result applies to classical equations of motion of the form ([9|) where the 
sum over y ranges over 7L d . 



2.6 Evolution of coherent states 



Denote by Si '■— sSi the unsealed spin operator in the spin-s-representation of su(2). For the polar 
angles (8, ip) £ [0, ir] x [0, 2tt) corresponding to the unit vector M £ E> 2 we define the coherent state in 
as 



C 2s+1 



\M) := cxp-^= 



\s) 
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where \s) is the highest-weight state, i. e. S z \s) = s \s). Note that \M) = e la ' s |s), where a = On and 
n is the unit vector (siny, — cosy, 0). 

e lv S- - e~ lip S+ 



Set A := ^ 



and U := e A so that \M) = U\s). Then using 

oo 1 

U*SiU = J2^l--[SuA],...,A], 

wc find 



k\ 

k=0 



U* SiU = sine cos ip S z + -^= cos 2 -(S++S-)--^ sin 2 - (er** S+ + eF* S-) , 

V 2 2 y2 2 

U*S 2 U = sm9sincpS z + cos 2 ^ (S+ - sin 2 °- {e 2lv S+ - e~ 2llp S^) , 

V2i 2 -y^i 2 

[TS 3 [/ = cos 9 S z - -= sine (eT*PS+ + e llp S-) . (19) 
v2 

As a consequence note that 

(sin cos yX 
sin(9siny J — M . (20) 
cos(9 / 

In order to derive our main result for coherent spins states we need the following lemma, which 
follows from direct calculations. 

Lemma 4. For any unit vector Mg§ 2 , we have that 

(M,S il ---S i M)-M il ---M ij 



Now let M : Z d — > § 2 be a configuration of classical spins on the lattice. Then M defines a state 
Pm on 21 as follows. For finite A, consider the product state 

\Ma) ■= (g) \M(x)) e J4?a ■ 

Then, for A € 21a, we set 

PM {A) := {Ma , A Ma) , 

and extend the definition of pm to arbitrary A £ 21 by continuity. 

Let M : R x Z d — > S 2 be the solution of the Hamiltonian equation of motion (|7|) with initial 
conditions M(0,x) — M(x). The following result links the quantum time evolution for coherent spin 
states with the corresponding classical configuration in the large-spin limit. 

Theorem 3. Let t E R, A e and M :Z d ^S 2 . Then 

lim pm^A) = A(M(t)) , 

s — >oo 

uniformly in t on compact time intervals. 

Proof. The proof is essentially a corollary of Lemma 0] and the proof of Theorem [TJ First, let \t\ < 
||V|| -1 . We know from $T3$ and (HU) that 

00 j-i 



1=0 

as well as 



A(M(t)) = £l lim {H A ,Af\M), 



pm&aA) = j2j\ ^ i p m ^ ^ (0 ) 

2=0 
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Now Lemma @] implies that 

\p M (M^)-M a \ < H^|. 
Arguing exactly as in the proof of Theorem [TJ we get the bound 

\ PM {a t A)-A{M(t))\ < (21) 

Arbitrary times are reached by iteration as in the proof of Theorem [TJ □ 

3 Mean-Field Limit 

This section is devoted to the dynamics of a quantum spin system in the mean-field/continuum limit. 
More precisely, we consider a system of quantum spins on a lattice with spacing h > 0. The limit h — > 
yields again a Egorov-type result: The quantum dynamics approaches the dynamics of a classical spin 
system defined on a continuum set. As in the previous section, we also discuss the thermodynamic 
limit and the time evolution of coherent states. 

3.1 A system of quantum spins on a lattice 

Let A C R d be bounded and open. We associate with each spacing h > the finite lattice 

A<*> := KL d n A. 

At each lattice site x € there is a spin of (fixed) magnitude s. The Hilbert space of this quantum 
system is 

Jf} h} := (g) C 2s+1 . 

The algebra of bounded operators on Jif^ is denoted by Sly^ . 

The spins are represented on Jtfg 1 ' by a family {Si(x) : i = 1, 2, 3, x G A^} of operators, where 
Si(x) is the i'th generator of the spin-s-representation of su(2), rescaled by h d /s. As usual, we replace 
the operators (Si,S 2 ,S 3 ) with (§ + ,S z ,S-). They satisfy the bounds ||5±|| < h d and \\S Z \\ = h d if 
s > 1, as well as \\S±\\ = \plh d and H-S^H = h d if s = 1/2. The commutation relations now read 

[Siix^Siiy)] = ^e ijk S(x,y)S k (x), (22) 

with i, j, k G /. 



3.2 A continuum theory of spins 

Let A C R d be a bounded, open set. A system of classical spins on A is represented in terms of the 
Poisson "phase space"H 

S A := {Me L°°(A;R 3 ) : \\M\U < l} , 

which we equip with the L°°-norm. In analogy to Section [21 we use the complex coordinates 
(M + , M z , M_) instead of (Mi, M2, M3), so that the Poisson bracket on Sa satisfies 

{M,(a;),M,(y)} = i e ijk S(x - y) M k (x) , (23) 

for i, j, fee/. 

3 As in the previous section, one may introduce a symplectic phase space F\ consisting of all M £ Sa such that 
\M(x)\ = 1 a.e. 
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In order to describe a useful class of observables on Ha, we introduce the space S§^ p \ p £ N, which 
consists of all functions / in C(W d ;C 3F ) that are symmetric in their arguments, in the sense that 
Pf = f, where 

(P f)ii...i p (xi, . . . , X p ) ■— — y ' /io.(i)... i„( p ) { x cr{l)i ■ ■ ■ ) x a(p)) ■ 



On the space we introduce the norms 



ll/H W : = hP d E \fn..,M,---^v)\> 

ii,...,i p el xi,...,x p GhZ d 

\\f\&x-=snp J2 hiP ' 1)d E \fi,...i P {x,x 2 ,...,x p 

il,---,ip€.I x 2 ,...,x p ehl, d 

We shall be interested in observables arising from / € SS^ satisfying 

l (ft) 



lim sup ||/||i < 00 ■ ( 24 ) 

h->0 

Note that Fatou's lemma implies that ||/||i < limsup ft _,. ||/||i • 

We define *Pa as the "polynomial" algebra of functions on 2a generated by functions of the form 

m a{I) '■= E / dx!---dx p fi 1 .. Ap (x 1 ,...,Xp)M il (x 1 )---M ip (x p ), 
— . JAp 

where / € 2%^ satisfies <Pa is clearly a Poisson algebra. We equip it with the norm H-AHoo := 

sup Mg ^ A |A(Af)| so that 

II^a(/)|U < ll/lli. (25) 

3.3 Quantization 

For / G @W let us define 

§A.(/) := E E fii-i p (xi,... > x p )S il (xi)---S ip (x p ). (26) 

ii,--')*j> x 1 ,...,x p £A( h 'i 

If / satisfies (|2"3|). we find that 

I|5a(/)|| < (27) 

As above, quantization ~ : <}3a — * 21 ^ is defined by M\(f) — : S\(f) : and linearity. Here : • : 
denotes the normal-ordering of the spin operators introduced above. Also, we set 1 = 1. Again, 
(A)* = % 

3.4 Dynamics in the mean-field limit 

We consider a family V = (V^™)) —1 of functions, where S satisfies 



V iCli n (xi,...,x n ) = y^] In {xi, . ■ ■ , x n ) , 
where, we recall, ~ on / maps (+, z, — ) to (— , z, +). We define the Hamilton function on Sa through 



H A := £ Af A (^ (n) ). (28) 

Set 



n=l 
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We impose the condition ||V|| := limsup^ ||V||^ < oo . In the continuum limit, we observe that 

oo „ 

Vne"sup V / dx 2 ---dx n \V^ ) An (x,x 2 ,...,x n )\ < \\V\\ , (30) 

, X.J 



as can be seen using Fatou's lemma. It now follows easily that, for each bounded set A, the sum ([28 
converges in || ■ on Sa and yields a well-defined real Hamilton function Ha- 
The Hamiltonian equation of motion reads 



d_ 

dt 



Mj(t,x) = i n J dx 2 ■ ■ ■ dx n V^l in (x, x 2 , 



n—l ii,...,i n ,j 



e tlij Mj (t, x) M l2 (t,x 2 )--- M in (t, x n ) . (31) 
By standard methods, we find the following global well-posedness result for (f3"Tj) . 

Lemma 5. Let A C K d by any open subset of M. d and Mo G 3a. Then (|31[) has a unique solution 
M G C (M, Sa) that satisfies M(0) = Mo- Moreover, the solution M depends continuously on the 
initial condition Mo- Finally, we have the pointwise conservation law \M(t,x)\ = \M(0,x)\ for all t. 

Remarks. 1. As in Section we denote the norm-preserving Hamiltonian flow by </> A . 
2. Time-dependent potentials V(t) may be treated exactly as in the previous section. 

Example. Consider 

H A = - dx h(t,x) ■ M(x) - - / dxdy J(x,y) M{x) ■ M(y) . 



Ax A 



which yields the Landau-Lifshitz equation of motion 
d 



-M(t,x) = M(t,x) A 



h(t,x)+ / dy J(x,y)M(t,y) 

J A 



Remark. In a formal way, the Landau-Lifshitz equation ^ mentioned in the introduction can be 
obtained from ([3]) with J — J(\x — y\) by Taylor expanding M (t, y) up to second order in y — x. This 
leads to ([2]) after rescaling time by 1 1— > at where a= ^ J J(|a;|)|a;| 2 dx. 

The quantum dynamics is generated by the Hamiltonian Ha £ 21^ defined as the quantization of 
H\. More precisely, each term of H\ is quantized and it may be easily verified that the resulting series 
converges in operator norm. The fact that H A is real immediately implies that H\ is self-adjoint. As 
above we introduce the short-hand notation 

a A A := Ao0 A , A e 21a, 

a h A := U h {t;H A )*AU h (t;H A ), A e ^ . 

Here, Uh(t; H A ) is the quantum mechanical propagator, equal to e lsh Hht if H A is time- independent. 

We are now in a position to state our main result on the mean-field dynamics of the quantum 
system on the finite lattice in the continuum limit, as h — * 0. 

Theorem 4. Let A C R d be open and bounded, A G ^Pa md e > 0. Then there exists a function 
A(t) G *Pa such that 

sup||a A A-^(i)||oo < e, (32) 

teR 

and, for any tgl, 

\\a\A- A(t)\\ < e + C{e,t,A)h d , (33) 
where C(e,t,A) is independent of A. 
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Proof. One finds, for / G 3§ {p) and g £ 3S^ , 

{M A (f),M A (g)}=pqM A (f^g) 
where / ge .5§ {p+q - 1) is defined by 



(34) 



(/ g)ix...i p+q -x{xi, ■ ■ -,Xp+q-l) 

■— iP ^ t Eijii fi'h-.-ip i&lt • • • x p) 9jip+i---i P +q~i (pi l 2-p+l j • • • j x p+q— l) ■ (35) 



We have the estimate 
where 



11/ -.g||i < ||/||=o4 Nli, 

,i := sup ^ I dx 2 ■ ■ ■dx p \f ll .. Ap (x, 

X , J 



(36) 



X2 1 • • • 3>T) 



Without loss of generality, we assume that A = M\(f) for some / € satisfying the bound 

Iterating 

00 

{H A ,M A (f)} = ^n P M A (y(™)^/) 

n=l 



we obtain that 

00 

{H A ,M A {f)} {l) = [pn x ][(p + n l -l)n 2 \---[(p + n 1 + ---+ni- x -l + l)ni] 

rti,...,raj=l 



M A (V (n!) (y^i-o (^< ni ) — /))), 



with norm 



{jr A ,M A (/)} (,) || oo < ^ [pni] [(p + ni- l)na] •••[(? + ni + .-. + m-i -I + l)m] 



"1 



< 



< e p ||/||in _ 
by Jill). Therefore, for |t| < ||V|| _1 , the series 



||^" ! )|U,i---||^ ( " l) ||cx )> i||/lll 

n £ (P + m + H ' ' + ni)l ni - m |U,i ■ • ■ nv^lUi H/lli 



(1) 



2=0 



(37) 



(38) 



converges in || • to a A A. 

The quantum case is dealt with in a similar fashion, with the additional complication caused by 
the ordering of the generators {Si(x)}. This does not trouble us, however, as an exact knowledge of 
the ordering is not required. It is easy to see that, for / and g as above, 



ish 



S A (f),S A (g) 



is equal, up to a reordering of the spin operators, to pqS A (f — g). Iterating this shows that 

(ish- d ) 1 [H A ,A\ (l) 
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is equal, up to a reordering of the spin operators, to 



^ [pni] [(P + n i - 1)^2] ■ ■ • [(p + ni H h - 2 + l)nj] 

ni,...,n;=l 

Consequently an estimate analogous to ([57)) yields, for s > 1, 

||(i S /r d )' [Ha,A\ (1) \\ < eP\\f\\[ h) n(\\V\\ (h) Y, 
which readily implies the bound 

\\jtU8h-*?[H A ,2\® <e* ll/Hi*) jt(\t\\\V\\W) 1 . 



(V 



(m) 



/) 



(39) 



If s = 1/2, the first line of (f3"T|) gets the additional factor \[2 



1=0 



«iH hni+p 



This may be dealt with by 



replacing the factor (p + ni + ■ • • + ni) 1 in the second line of (f3"T|) with (rp + rn\ + • • • + rn{) 1 jr l . The 
desired bound then follows for < r < 1 — h log 2. Note that in this case the convergence radius for t 
is reduced to r||F|| _1 . For ease of notation, we restrict the following analysis to the case s > 1, while 
bearing in mind that the extension to s — 1/2 follows by using the above rescaling trick. 

Now, by definition of ||V||, for any \t\ < there is an ho such that (|39[) converges in norm to 

ofj^A for all h < Iiq, uniformly in h and A. 

In order to establish the statement of the theorem for short times \t\ < ||V|| _1 , we remark that 
the commutation relations (|22|) imply the bound 

\\A-B\\ < ^P 2 \\f\\?\ 

for arbitrary reorderings, A and B, of the same operator Sa(/), with / 6 3§&> for some p < 00. 
If we define a^A through its norm-convergent power series, we therefore get 

||3 A A- a A A|| 
h d ^ | t |J 



< — n E M ^ p + "1 - !) n 2] ■■■[{p + n 1 + --- + n l _ 1 - I + l)nj] 

1=0 ' m,...,ni 

(ni + ... + ni _ Z + l )2 ||y(n.)||W 11/11 W 

<-±\t\ i e (p+ " l+ ;r +n,),+a "i-"«ii^iigi- 



1^)11^11/11^ 



i=0 m,...,raj 



< -Ei i i ieP ii/iii' l) ( / + 2 )( ? + 1 ) E™ e "i^ 



n IIT/C" 1 )!!^) 

loo.l 



2=0 



< 



e p ii/iirE( ? + 2 )( ? + 1 ) (nii^ii 



(h)\J 



2=0 



where in the last step we have used the fact that the sum convergences uniformly in h, for h small 
enough, as seen above. 

Arbitrary times are reached by iteration of the above result. □ 
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3.5 The thermodynamic limit 

The above result may again be formulated in the thermodynamic limit as A — > KL d . We only sketch 
the arguments, which are almost identical to those of Section [2751 
The quantum quasi-local algebra is 

m := v 

AccR d 

The existence of dynamics is guaranteed by the following statement. 

Lemma 6. Let h > and suppose A G 21^ for some bounded and open Aq C R d . Then, for any 
t G R, the following limit exists in the norm sense: 

lim a\A =■ a* A, 

A^oo 

By continuity this extends to a strongly continuous one-parameter group (a')tgR of automorphisms of 
The classical quasi-local algebra is 

21 := \/ <Pa- 

ACCR d 

Lemma 7. Lei A G *Pa f or some open and bounded Aq CC M d . Then, for any t G M. i/ie following 
limit exists in \\ ■ W^: 

lim a A ;4 =: a*A, 

A^oo 

Zfy continuity this extends to a strongly continuous one-parameter group (a )tgR o/ automorphisms of 
21. Furthermore, a 1 A = A o 0', where $ = (j)^ d is the Landau- Lifschitz flow defined in Lemma\^ 

Now, for / G M(/) and are well-defined in the obvious way. Define as the algebra 

generated by functions of the form M(f), where / satisfies (|24|) . 

Theorem 5. Let A G and e > 0. XTien i/iere exists a function A(t) G swc/i i/iai 

BupHa^-ACtJHoo < e, (40) 

and, for any igl, 

||a'A-A(t)|| < e + C(e,t,A)h d . (41) 

3.6 Evolution of coherent states 

In this section, our "smearing functions" / are assumed to have compact support, i. e. to belong to 
the space 

#w : = ^'nCetr^c 3 '). 

In addition, we require the interaction potential V to be of finite range in the sense that there exists 
a sequence R n > such that if \xi — Xj\ > R n for some pair then V^ A (xi, . . . , x n ) = 0. 

Next, we take some initial classical spin configuration M G C(Mr; § 2 ), or, more generally, a function 
M : M. d — > § 2 whose points of discontinuity form a null set. We shall study the time evolution of 
product states pu on 21^ that reproduce the given classical state M. For open and bounded A C R rf , 
we define the product state 

|Ma):= (g) |M(x)>, 

where \M(x)} is the coherent spin state corresponding to the unit vector M(x). For A G 2t A , define 

Pm{A) := (M A , A M a ), 
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which we extend to arbitrary A S by continuity. 

For our main result on the time evolution of coherent states, we first record the following auxiliary 
result whose elementary proof we omit. 

Lemma 8. Let f £ SS^ satisfy (|24]). Then 

lim p M (£(/)) = M(f). (42) 

The last result in this paper links the quantum time evolution of coherent spin states with the 
classical evolution in the mean-field/continuum limit when the lattice spacing h tends to 0. 

Theorem 6. Let t £ R, A e *p and M be as described above. Let M(t) be the solution of (f5Tj) on M. d 
with initial configuration M . Then 

lim p M (a* A) = A(M(t)) , 
uniformly in t on compact time intervals. 

Proof. The proof is a corollary of the proof of Theorem [H First, let \t\ < and pick an e > 0. 

Choose a cutoff such that the tails of the thermodynamic limits of the series ([55]) and (j3"9")l are bounded 
by e. We therefore have to estimate a finite sum of terms of the form 

\p M (S(g))-M(g)\, 

where g S 3§i P< " 9 ^ because of our assumptions on V. By Lemma [SJ for h small enough, these are all 
bounded by e, and the claim for small times follows. Finally, by iteration, we extend the result to 
arbitrary times. □ 

Acknowledgement. Wc thank Simon Schwarz for very useful discussions. 
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